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ABSTRACT 


Recent observations of the high-mass X-ray binary Cygnus X-1 have shown that both the companion star (41 
solar masses) and the black hole (21 solar masses) are more massive than previously estimated. Furthermore, 
the black hole appears to be nearly maximally spinning. Here we present a possible formation channel for the 
Cygnus X-1 system that matches the observed system properties. In this formation channel, we find that the 
orbital parameters of Cygnus X-1, combined with the observed metallicity of the companion, imply a significant 


reduction in mass loss through winds relative to commonly used prescriptions for stripped stars. 
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1. INTRODUCTION 


Cygnus X-1 is a high-mass X-ray binary (HMXB) in the 
Cygnus OB3 association, which hosts a star in orbit with 
a black hole (BH) (e.g. Webster & Murdin 1972; Bolton 
1972, 1975; Hutchings et al. 1973). The BH accretes mat- 
ter from the stellar wind; this accretion powers X-ray radia- 
tion (e.g. Davidson & Ostriker 1973; van den Heuvel 1975; 
Conti 1978; Petterson 1978) and a jet (e.g. Bisiacchi et al. 
1974; Marti et al. 1996; Stirling et al. 2001). Orosz et al. 
(2011) inferred the BH and stellar companion masses of 
Cygnus X-1 to be 14.8 + 1.0Mọ and 19.2 + 1.9Mo, re- 
spectively. Revised measurements of the distance to Cygnus 
X-1 (Miller-Jones et al. 2021) indicate that both objects are 
significantly more massive. The temperature and luminos- 
ity of the optical companion are estimated to be Tef = 
31.1 + 0.7 kK and log(L/Lo) = 5.63 + 0.07 with a mass 
of Mopt = 40.6* 77 Mo, where we quote the median value 
and the 68 per cent confidence interval boundaries (Miller- 
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Jones et al. 2021). The mass of the BH is estimated as 
Mpg = 21.9722 Mo. The binary has an almost circular 
orbit with a semi-major axis a = 0:94 onis AU and ec- 
centricity e = 0.0189* 0.0075 (Miller-Jones et al. 2021). The 
BH is inferred to be nearly maximally spinning with a di- 
mensionless spin of at least 0.95 according to both disk con- 
tinuum and reflection line fitting studies (Gou et al. 2011; 
Fabian et al. 2012; Miller-Jones et al. 2021; Zhao et al. 2021). 
Via optical spectroscopy, it has been found that the ratios of 
the surface abundances of both helium and iron to hydrogen 
are about twice the respective values for the Sun (Shimanskii 
et al. 2012). As we discuss below, these observations taken 
together present a challenge for models of massive stellar bi- 
nary evolution. 

In this paper we describe the constraints that these obser- 
vations place on the Cygnus X-1 formation channel (Sec. 2). 
We explore how the helium main sequence (HeMS) phase 
of this channel provides a constraint on the wind mass-loss 
rates of massive stars (Sec. 3). We present the likely future of 
the system and its implications for gravitational-wave detec- 
tions (Sec. 4). Finally we discuss some caveats and questions 
raised by the observations of Cygnus X-1 (Sec. 5). 
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Figure 1. Assumed formation channel for Cygnus X-1. 1: Two 
stars (primary top, secondary bottom) are born in a binary. 2: The 
primary star is more massive and evolves faster, expands and starts 
stable mass transfer onto the secondary. 3: The primary is left as 
a hot stripped HeMS star with a companion which could have ac- 
creted a significant amount. 4: The primary star collapses and leaves 
behind a BH orbiting the secondary MS star. 


2. CYGNUS X-1: OBSERVATIONS AND ASSUMED 
CHANNEL 


We assume the following formation channel for Cygnus 
X-1 (see Fig.1 for illustration). Two stars are born in a bi- 
nary. The more massive star (the primary) evolves more 
quickly and expands first. The companion (the secondary) 
is close enough for the late main sequence primary to com- 
mence mass transfer. The primary is stripped of its envelope, 
leaving an exposed helium core. The core continues nuclear 
fusion until it collapses and forms a BH in orbit with the still 
core-hydrogen-burning main sequence (MS) companion. In 
the sections below, we describe the observations and theoret- 
ical analyses that support this channel. 


2.1. Eccentricity, peculiar velocity and fallback 


The collapse of a star into a compact object can impart both 
a natal kick from an asymmetric explosion (see e.g. Fryer 
2004 and references therein) and a kick from rapid symmetric 
mass loss to the system (Blaauw 1961). The low eccentricity 
of the Cygnus X-1 binary (Orosz et al. 2011; Miller-Jones 
et al. 2021) seems to disfavour a significant natal kick. How- 
ever, tidal forces could have circularised the system since the 
collapse: for the inferred stellar and binary properties, the 
circularisation timescale for dynamical tides (applicable for 
radiative-envelope stars, Zahn 1977) is only ~ 10? years (es- 
timated using equations (41) and (43) of Hurley et al. 2002). 

On the other hand, the system's small peculiar veloc- 
ity strongly indicates that the BH experienced nearly com- 
plete collapse with little mass ejection during its forma- 
tion (Mirabel & Rodrigues 2003). The peculiar velocity is 
10.7 + 2.7 km s-! relative to its host association (Rao et al. 
2019), which limits the amount of instantaneous symmetric 
mass loss (Blaauw 1961; Nelemans et al. 1999; Wong et al. 


2012) to S 2Mo, unless the kick from symmetric mass ejec- 
tion is fortuitously cancelled by an oppositely directed natal 
kick. Further indirect evidence for the low natal kick lies in 
the absence of Type C quasi periodic oscillations, which may 
indicate Lense-Thirring precession due to spin-orbital mis- 
alignment as a result of a natal kick (Stella & Vietri 1998) 
and are observed in most known black hole X-ray binaries 
(Ingram et al. 2016), the majority of which are believed to 
have received strong (~ 100 km s^ !) natal kicks (Atri et al. 
2019). The low amount of ejected mass is also consistent 
with the current eccentricity in the absence of significant 
tidal circularisation, and suggests nearly complete fallback 
onto the black hole except for a small amount of neutrino 
mass loss (e.g. Nadezhin 1980; Lovegrove & Woosley 2013; 
Fernández et al. 2018). Nearly complete collapse matches 
the theoretical models simulating the fall-back onto black 
holes of similar masses (Fryer et al. 2012) and observational 
evidence for massive stars disappearing without supernovae 
(Adams et al. 2017). 


2.2. Black-hole progenitor mass 


The Eddington limit on the mass accretion rate for a BH 
of this mass is < 2 x 107" Mc yr-!. Under the assump- 
tion of Eddington-limited accretion (but see, e.g., Fragos & 
McClintock 2015; Eldridge et al. 2017; van Son et al. 2020, 
who relax this assumption), only a negligible amount of mass 
could have been accreted onto the BH since it formed (King 
& Kolb 1999). The companion lifetime sets the accretion 
duration to no more than a few Myr, meaning that at most 
< few x 107! Mc could have been accreted. Assuming that 
the jet turned on promptly after the formation of the black 
hole, the few x 10^ yr estimate of the age of the jet (Russell 
et al. 2007) places an even stronger constraint on the amount 
of accreted mass, < 107? Mo (see also Sell et al. 2015, who 
consider other models for the nebula origin but reach simi- 
lar conclusions about its age). Therefore, the current black 
hole mass Mgy = 21.2722 Mọ is a good estimate for the 
progenitor mass just before the collapse. 


2.3. Black hole spin & progenitor spin 


The BH has a dimensionless spin x close to unity, x > 
0.95 (Gou et al. 2011; Fabian et al. 2012; Miller-Jones et al. 
2021), although spin measurements may have large sys- 
tematic errors because of modelling assumptions (Miller & 
Miller 2015; Kawano et al. 2017; Zhao et al. 2021). A spin 
close to unity, coupled with negligible mass loss during col- 
lapse, suggests that if the pre-collapse progenitor had excess 
angular momentum, x > 1, it must have been carried away 
by a small amount of mass with high specific angular mo- 
mentum (cf. Janka 2013, 2017; Batta & Ramirez-Ruiz 2019; 
Murguia-Berthier et al. 2020). 

The rotational angular momentum of the BH was either 
present in the progenitor, or was gained during or after its 
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collapse. Here we briefly summarise why we assume that the 
angular momentum comes from the progenitor and why this 
implies that the progenitor must have been stripped early in 
its evolution (see Mandel & Fragos 2020 for a longer discus- 
sion). 

A BH needs to roughly double its mass through accretion 
in order to go from zero to maximal spin (Bardeen et al. 1973; 
Thorne 1974). This is not possible for Cygnus X-1 under the 
assumption of Eddington-limited accretion. The spin could 
have been acquired during the collapse, e.g., if the compan- 
ion torques some of the ejecta which then fall back onto the 
BH (Batta et al. 2017; Schrøder et al. 2018), but this requires 
some fine-tuning: for example, the ejecta must have suffi- 
cient velocities to be torqued by the companion, but with- 
out a significant tail of escaping ejecta to match formation 
through nearly complete fallback as discussed above. There- 
fore we assume that the angular momentum was present in 
the progenitor at the moment of the collapse. 

Observations suggest that stars might be born as rapid ro- 
tators (Fukuda 1982; Rosen et al. 2012; Ramirez-Agudelo 
et al. 2013, 2015), although these observations could be af- 
fected by binary interactions (Langer et al. 2008; de Mink 
et al. 2013). Even if stars are born as rapid rotators, they 
spin down through wind-driven mass loss and are unlikely to 
retain enough angular momentum to form rapidly spinning 
BHs. 

The black hole’s dimensionless spin is determined by the 
ratio of its total angular momentum to the square of its mass. 
At lowest order, it is therefore insensitive to the redistribution 
of angular momentum through the progenitor star. However, 
if the bulk of the angular momentum moves into the enve- 
lope as the star expands, and the envelope is subsequently 
removed by winds or mass transfer, the remaining core is 
likely left with too little angular momentum to produce a 
rapidly spinning BH (Petrovic et al. 2005; Belczynski et al. 
2017; Fuller & Ma 2019; Bavera et al. 2020; Mandel & Fra- 
gos 2020). On the other hand, the core could retain sufficient 
angular momentum to form a rapidly spinning BH if there is a 
large amount of differential rotation between the layers of the 
star (Hirschi et al. 2005). Therefore, the efficiency of angular 
momentum transport in the star plays a key role, yet both the 
mechanism and degree of core-envelope coupling remain un- 
certain. Theory (Tayler 1973; Spruit 2002; Fuller et al. 2019; 
Takahashi & Langer 2020) and observations, such as rotation 
rates of low-mass giant stars (Cantiello et al. 2014), suggest 
that there may be efficient angular momentum transport and 
the envelope is coupled to the core. If so, the observed rapid 
spin of the Cygnus X-1 black hole in a close binary, in which 
the progenitor must have lost its envelope, appears to require 
pre-collapse interaction with the companion to spin up the 
stellar core. 


Tidal locking of the period of the stellar rotation to the 
period of the binary provides the most likely mechanism 
for producing a rapidly rotating BH progenitor (Izzard et al. 
2004; Kushnir et al. 2016; Zaldarriaga et al. 2018; Belczyn- 
ski et al. 2017; Bavera et al. 2020). Chemically homogeneous 
evolution could yield rapidly rotating black holes (Mandel & 
de Mink 2016; Marchant et al. 2016), but is not expected to 
operate at such high metallicities and is not consistent with 
the observed expansion of the companion. 

Instead, the following evolutionary sequence, proposed by 
Valsecchi et al. (2010) for M33 X-7, a similar high-mass X- 
ray binary with a rapidly spinning BH, and investigated in 
detail by Qin et al. (2019), appears to be the most likely for- 
mation mechanism for Cygnus X-1. The binary starts out 
with a period somewhat shorter than the current observed 
one. The more massive primary — the BH progenitor — com- 
mences mass transfer while still on the main sequence. This 
mass transfer, is likely to be largely non-conservative if lim- 
ited by the spin-up of the accretor (Qin et al. 2019); this is 
favoured both by the observed evolutionary state of the sec- 
ondary (conservative mass transfer would imply an initially 
lower-mass primary, and hence an older binary, placing an 
upper limit on the BH progenitor mass that would make it un- 
likely to form such a massive BH) and the observed period. 
The mass transfer removes the donor's envelope, prevent- 
ing subsequent re-expansion and angular momentum loss. 
Meanwhile, the core remains tidally locked on the MS. Af- 
ter hydrogen is exhausted in the core at the end of the MS, 
the star contracts into a rapidly spinning HeMS star. While 
this HeMS star is no longer tidally locked, it can still col- 
lapse into a rapidly spinning BH. Qin et al. (2019) find that 
the efficiency of angular momentum transport does not play a 
significant role for the evolution of the black hole progenitor 
during the main sequence, where tidal locking keeps its core 
rapidly spinning, but could be a key factor in determining the 
ultimate black hole spin through its impact on the angular 
momentum lost through winds in later evolutionary stages. 
The Qin et al. (2019) models successfully reproduce systems 
with the orbital parameters of Cygnus X-1. In fact, the match 
appreciably improves with the upward revision in the BH and 
optical companion masses, as the observed masses in the bot- 
tom left panel of figure 3 of Qin et al. (2019) shift toward the 
locus of their model evolutionary trajectories (though their 
evolutionary trajectories are only shown for an initial mass 
ratio of 0.4, while the latest observations support more com- 
parable masses). Given the number of uncertainties relating 
to this first mass transfer episode, in what follows we focus 
on the subsequent evolution of the binary under the assump- 
tion of this formation channel for Cygnus X-1. 


2.4. Companion mass & abundance 
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We compared the observed mass, luminosity and temper- 
ature of the optical companion against analytic fits to stellar 
tracks of Hurley et al. (2000) as implemented in the COM- 
PAS rapid population synthesis code (Stevenson et al. 2017; 
Vigna-Gómez et al. 2018). The observations are consistent 
with a MS star that is about 70 to 80 per cent through its core 
hydrogen burning phase, ignoring the impact of rotation. 

These stellar tracks are for stars with regular hydrogen-rich 
atmospheres. Using them ignores the possible impact of non- 
standard surface abundances, and so corresponds to the as- 
sumption that only a thin surface layer has a significant over- 
abundance of helium, rather than a uniform distribution of 
enriched material throughout the star. Stars in later stages of 
the MS with the relevant mass and metallicity should have at 
most a very thin convective layer at the surface (e.g., Maeder 
et al. 2008; Kippenhahn et al. 2012). Therefore, mixing is 
expected to be relatively inefficient: in the absence of large- 
scale convection, the Rayleigh-Taylor instability is likely to 
be suppressed by temperature inversion in the accreted ma- 
terial (Kippenhahn et al. 1980; Braun & Langer 1995). The 
resulting thermohaline mixing will mix the helium-rich ma- 
terial through the companion only on timescales longer than 
the expected few x 10* years since the formation of the BH if 
the bulk of the enriched material was accreted at or shortly 
before the BH formation. 

The evolutionary channel shown in figure 1 assumes that 
the optical companion has not overflowed its Roche lobe. 
While previous studies explored this possibility, perhaps with 
intermittent Roche lobe overflow followed by longer periods 
when the binary is detached as in the present state, these 
models were typically based on assumptions that the com- 
panion is less massive than the BH (Podsiadlowski et al. 
2003), which is inconsistent with present observations. In- 
deed, if we assume that mass transfer onto a black hole is 
almost entirely non-conservative because of the Eddington 
limit, the binary's semi-major axis a evolves as 


a = 9 Mot 1 E 1 Mopt 
a Mox MET Mox + Mpu 


| , Q.D 


where y is the specific angular momentum of the ejected ma- 
terial in units of the binary's specific orbital angular momen- 
tum J/(Mopt + Afen), The binary can widen as a result of 
such mass transfer only if y < 1 for the observed component 
masses. However, in the common assumption of isotropic 
re-emission from the accretor, y = Moy /Mpu ~ 2. The 
ejected material would have to carry much lower specific an- 
gular momentum in order for the binary to be able to detach 
once Roche lobe overflow from the companion commences, 
which seems unlikely, lending support to our proposed chan- 
nel. 

The estimated mass and lifetime of the companion already 
enable us to roughly infer the initial mass of the BH progen- 


itor. We assume that both stars are born and start fusion at 
the same time and use a fit (Farr & Mandel 2018) to the Brott 
et al. (2011); Kóhler et al. (2015) stellar models for the MS 
lifetime of non-rotating massive stars. For this simple esti- 
mate, we ignore the impact of binary interactions, consistent 
with the assumption of largely non-conservative mass trans- 
fer. The BH progenitor should have had an initial mass of 
~ 55 — T9 Mg, in order to complete its evolution while leav- 
ing behind a MS companion of mass Mp4 at ~ 70% — 80% 
of its MS lifetime. 

The surface abundances of the companion are non- 
standard for massive MS stars and a challenge to explain even 
in the context of binary interaction. We discuss the abun- 
dances in more detail in Section 5. Here, we focus on the 
iron abundance, as this is key to the analysis of line-driven 
winds in the following section. Shimanskii et al. (2012) 
find that the iron abundance of the companion is 2.2 times 
the solar iron abundance, although precise measurements are 
challenging due to the complexity of the system. Assum- 
ing Zo = 0.014 for solar metallicity (Asplund et al. 2009), 
this corresponds to an effective metallicity of Z ~ 0.03. On 
the other hand, Daflon et al. (2001) find a slightly sub-solar 
iron abundance of log e(Fe) = 7.33 + 0.12 in HD 227460, 
which is a BO.5V star in the same Cygnus OB3 association 
as Cygnus X-1. This could indicate that Z ~ 0.01 is a better 
estimate of the initial iron abundance, and the iron abundance 
of the companion in Cygnus X-1 has been enhanced during 
the collapse of the primary to a BH. Therefore, we explore 
an initial metallicity range 0.01 < Z < 0.03 in the following 
section. 


3. CYGNUS X-1: MAXIMUM WIND MASS-LOSS RATE 


Hereafter we assume that, after the mass transfer episode 
induced by the BH progenitor, the primary is left with no 
hydrogen layer on top of the He core. This assumption is 
consistent with the channel proposed by Qin et al. (2019), in 
which the rapidly spinning BH is ultimately formed through 
the collapse of a Wolf-Rayet star (see however discussion in 
Section 5). 

In order for Cygnus X-1 to form through the channel de- 
picted in Fig. 1, two things must be true. Firstly, the HeMS 
star must be born with sufficient mass to give rise to the ob- 
served BH mass even after losing mass through Wolf-Rayet 
winds during the HeMS. However, at a given metallicity, 
there is a maximum to the He core mass that can be formed 
at the end of the MS: as more massive stars have higher wind 
mass loss rates, terminal-age MS He core masses asymptot- 
ically approach a maximum as a function of zero-age MS. 
This maximum He core mass is plotted in figure 2. Be- 
cause the BH progenitor is stripped of its hydrogen enve- 
lope at the end of the MS for the binary evolutionary channel 
we assume, we can directly constrain the amount of mass 
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Figure 2. Maximum helium core mass at terminal-age MS as a 
function of metallicity, maximized over zero-age MS masses, based 
on the stellar tracks of Hurley et al. (2000) as implemented in COM- 
PAS (Stevenson et al. 2017; Vigna-Gómez et al. 2018). 


loss during the HeMS phase. The Wolf-Rayet winds must 
then not remove more than the difference between this max- 
imum mass and the final BH mass. We refer to this as the 
“MyHems € MneMsS,max " condition. 

The maximum He core mass plotted in figure 2 is very sen- 
sitive to the MS wind prescriptions (see discussion in Renzo 
et al. 2017; Neijssel et al. 2019; Miller-Jones et al. 2021). The 
analysis of the remnant mass from massive single stars is par- 
ticularly sensitive to luminous blue variable winds (Belczyn- 
ski et al. 2010; Miller-Jones et al. 2021), which can remove 
the hydrogen envelope of the star and hasten the onset of the 
Wolf-Rayet phase (Conti 1975); therefore, increasing the po- 
tential remnant mass can be achieved by decreasing either 
luminous blue variable winds or Wolf-Rayet winds. The plot 
in figure 2 assumes the default COMPAS luminous blue vari- 
able mass loss rate of 1.5 x 10~4Mo yr-! (Belczynski et al. 
2010) for stars approaching the Humphreys-Davidson limit 
(Humphreys & Davidson 1994). The convective overshoot- 
ing parameter employed in stellar evolution calculations is 
another source of uncertainty, with greater overshooting lead- 
ing to larger cores (Brott et al. 2011). To explore the im- 
pact of convective overshooting, we simulated the He core 
mass at terminal-age MS for stars with different metallicities 
starting from a zero-age MS mass of 150 M, using the stel- 
lar evolution code MESA (v12115, Paxton et al. 2011). We 
find that varying the overshooting parameter by an order of 
magnitude between f — 0.02 and f — 0.2 (Qin et al. 2019 
used f — 0.11) with the step overshoot scheme changes the 
He core mass at terminal-age MS by between 196 and 1396 
for the range of metallicities in figure 2. Finally, the COM- 
PAS single stellar evolution models are based on the fitting 
formulae of Hurley et al. (2000) to the evolutionary tracks 
of Pols et al. (1998) which involve extrapolations to higher 
masses than those covered by the initial range of models, and 
may under-estimate the He core masses of massive stars. We 
therefore use the following additional constraint, which by- 
passes these sources of uncertainty in the maximum He core 
mass. 


In our assumed channel the secondary MS companion 
must not overflow its Roche lobe onto the HeMS primary. 
We write this constraint as (Eggleton 1983) 


0.49q?/3(t) 


Ra(t) € a(t) 0.642/3 (t) + In(1 + q'/3(t))' 


(3.1) 


where mə2(t) and R2(t) are the mass and radius of the com- 
panion as a function of time, m; (t) the mass of the BH or 
its progenitor, q(t) = IL is the mass ratio, and a(t) the 
orbital separation of the binary. 

The orbital separation widens due to wind mass loss. In the 
limit of fast, non-interacting winds, the widening is described 


by 


à ` Ma 

——— 32 

" Mas (3.2) 
where Mio = m4 + ma is the total mass of the system. 


Because winds widen the binary and remove mass from the 
HeMS primary faster than from its MS companion, they in- 
crease the size of the secondary's Roche lobe over time. 
Thus, even though the secondary is not overflowing its Roche 
lobe now, it may have done so in the past if the mass-loss 
rate was high. If we evolve the system back in time, the re- 
quirement that the secondary never overflows its Roche lobe 
imposes an alternative upper limit (Axelsson et al. 2011) 
on the maximum Wolf-Rayet wind mass-loss rate. We re- 
fer to this as the "no Roche-lobe overflow (RLOF)" condi- 
tion. Although the non-interacting wind assumption, describ- 
ing the widening of the binary (Equation 3.2), may be an 
over-simplification for such short-period systems (MacLeod 
& Loeb 2020), it allows us to conservatively estimate the 
constraints imposed by the existence of Cygnus X-1 and pre- 
sented below. 

We parametrise the mass-loss rate through Wolf-Rayet 
winds, modelled with the prescription proposed by Belczyn- 
ski et al. (2010) and based on (Hamann et al. 1995; Hamann 
& Koesterke 1998; Belczynski et al. 2010), with a multiplica- 
tive parameter fwr, following Barrett et al. (2018) (see Ap- 
pendix A for a definition and discussion). We constrain the 
allowed parameter space of the wind strength by rewinding 
the evolution of the binary from the current state. Because the 
BH formed recently in our model, we set the luminosity and 
temperature of the MS companion at the end of the HeMS 
phase of the primary equal to the current inferred luminos- 
ity and temperature of the observed MS secondary. The mass 
and age of the secondary inferred from temperature and lumi- 
nosity vary slightly for different metallicities. As we argued 
earlier, the BH is expected to lose negligible mass during col- 
lapse, so we set the mass of the HeMS primary at the end of 
that phase equal to the inferred BH mass. We assume that 


6 NEIJSSEL ET AL. 


the secondary was 99.7% Roche-lobe filling at the end of the 
primary's HeMS phase (Miller-Jones et al. 2021)! . 

The reverse evolution of the MS and HeMS stars is fol- 
lowed using the analytic fits to the stellar tracks of Pols et al. 
(1998) as presented in Hurley et al. (2000). The winds of the 
MS star are given by Vink et al. (2001). The HeMS winds 
are parametrised with the multiplicative factor fwg as de- 
scribed above and in Appendix A. The orbital response to 
mass loss is given by Eq. (3.2). We go back in the evolu- 
tion for a HeMS lifetime (note that the HeMS lifetime de- 
pends on how massive the HeMS star initially was, which 
depends on the wind strength we assume) and check that 
the Myems < Muems,max condition is satisfied at the start 
of the HeMS phase and the no RLOF condition is satisfied 
throughout this phase. 

Figure 3 shows the upper limits on the wind strength, 
parametrized as fwr, imposed by these two conditions as 
a function of metallicity. Both conditions show that the wind 
strength has to be reduced from the nominal value fwr = 1 
throughout the range of metallicities we have explored. The 
strongest constraint is placed by the Myems < Muyems,max 
condition, which is subject to uncertainties in the MS wind 
strengths, overshooting, and the single stellar evolution fits of 
Hurley et al. (2000). This yields an upper limit fwr < 0.4 
at Z = 0.01 and < 0.05 at Z = 0.03. However, even if 
we lift this constraint, the no RLOF condition still places a 
strong constraint on the allowed mass loss rate: fwr < 0.45 
at Z = 0.01 and S 0.15 at Z = 0.03. We also explore 
the impact of BH mass and Roche-lobe filling factor mea- 
surement uncertainties and find that the constraints on fwr 
change by < 0.1 over the range consistent with observations 
(Miller-Jones et al. 2021). 


4. CYGNUS X-1: FUTURE EVOLUTION 


The revised mass of the companion star in Cygnus X-1 
makes it a potential candidate for a future BH and raises 
the intriguing prospect that this system could form a merg- 
ing binary BH, connecting HMXBs with gravitational-wave 
sources. We model the future evolution of the system and 
find that it is unlikely to form a merging binary BH. 

Belczynski et al. (2011) used earlier, lower estimates of 
the mass of the BH and companion in Cygnus X-1 in order 
to analyse the future evolution of this system. They predicted 


! In Miller-Jones et al. (2021) the Roche-lobe filling factor is defined as the 
ratio between the distance from the center of the star to the point where 
the equipotential surface enclosing the volume of the star crosses the axis 
connecting the star to its BH companion, and the distance from the center of 
the star to the L1 Lagrange point. In this paper, we define the Roche-lobe 
filling factor as the ratio between the radius of the star and the volume- 
equivalent Roche-Lobe radius of Eggleton (1983). We therefore convert 
the median and one-c lower bounds of 0.96 and 0.93 reported in table 1 of 
Miller-Jones et al. (2021) to values of 0.997 and 0.99, respectively, that we 
use here. 


no RLOF 
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Figure 3. Upper limits on the parametrised Wolf-Rayet wind mass 
loss rate from HeMS stars. The Myems < Muems,max (purple) 
and no RLOF (teal) conditions only allow fwr values below the 
curves at metallicity Z. The dashed curves correspond to assuming 
An = 19.2Mo and the current companion Roche-filling factor 
of 0.990 rather than the default values of 21.2 Mo and 0.997 (solid 
curves), and indicate the impact of observational uncertainty. 


that imminent mass transfer from the nearly Roche-lobe fill- 
ing companion will significantly reduce the companion mass, 
leaving behind a core that could only form a neutron star, 
not a BH. They further estimated that the natal kick has a 
30% probability of unbinding the binary, and there was only 
a ~ 1% probability that the ensuing neutron star — black hole 
binary would merge within 14 Gyr through gravitational- 
wave emission. 

Here we update their predictions based on revised observa- 
tions, using the COMPAS population synthesis code (Steven- 
son et al. 2017; Vigna-Gómez et al. 2018) for binary evolu- 
tion calculations. We estimate the initial masses and orbital 
separation as in Section 3. We assume a Roche lobe filling 
factor of 0.997 and a metallicity of Z — 0.02. The median 
values for the BH mass and the luminosity and temperature 
of the MS companion from Miller-Jones et al. (2021) then 
yield Mpy = 21.2Mo, Mopt = 38.9Mo, a = 51.3 R5, and 
TMs = 0.81, i.e., the MS companion is 81 per cent through its 
core hydrogen burning phase by time. Once the MS compan- 
ion overflows its Roche lobe, the mass transfer phase could 
brighten Cygnus X-1 up to close to its Eddington luminos- 
ity of a few x 10?? erg s^!. Once the envelope is removed, 
the companion will appear as a Wolf-Rayet (WR) star in an 
HMXB in a phase lasting for ~ 3 x 10° years. In the de- 
fault model described below, we apply WR mass loss with 
fwR = 1 to such a star. 
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If we follow a treatment similar to Belczynski et al. (2011), 
in which the He core mass is determined by the mass of 
the star at the end of MS (after stripping in this case), we 
find the companion will collapse into a neutron star. Even 
if the binary survives the neutron star natal kick, which hap- 
pens in 7% of all binaries in our models, it will generally 
be too wide to merge through gravitational-wave emission, 
with only a 0.12% probability that it will merge in 14 Gyr, 
similar to the findings of Belczynski et al. (2011). During 
this time, it may be detectable in radio pulsar surveys as a 
neutron-star — black-hole binary, although the non-recycled 
pulsar will likely only be observable as a radio source for a 
few tens of Myr. However, this treatment may significantly 
under-predict the core masses of stars that donated mass on 
the MS. 

In our revised model, we address the potentially under- 
estimated core mass of stripped MS donors under the as- 
sumption that this core mass is determined only by the mass 
of the star at the end of the MS. We account for the substan- 
tial amount of helium synthesized by the companion before 
interaction with the following crude approximation. We con- 
sider a constant rate of helium production during the MS, 
so the mass of helium produced is Myems = TusMe,final, 
where Me,final is the final helium core mass that would be 
achieved at the end of the MS phase for the given stellar mass, 
in the absence of mass transfer. Keeping track of the helium 
mass synthesized before stripping, Min, us, leads to higher 
core masses at the end of the MS phase. 

With this correction, we predict that in the absence of BH 
natal kicks, Cygnus X-1 will become a bound binary BH, 
but one that is too wide to merge within 14 Gyr. If we in- 
corporate the COMPAS prescription for BH kicks based on 
the "delayed" model of Fryer et al. (2012), this changes to a 
38% probability of surviving the kick and forming a bound 
binary BH, with a 496 probability that this binary will merge 
within 14 Gyr through gravitational-wave emission due to a 
fortuitous kick. However, the kick prescription for low-mass 
black holes that do not undergo complete fallback (Fryer 
et al. 2012) is rather uncertain, with conflicting evidence on 
the natal kick magnitudes of low-mass BHs (Repetto et al. 
2017; Mandel 2016; Mirabel 2017; Wyrzykowski & Mandel 
2020; Atri et al. 2019). 

We also consider the impact of varying WR mass loss. Mo- 
tivated by the results reported in figure 3 for Z — 0.02, we 
reduce fwr to 0.2 from the value of 1 considered above. This 
increases the remnant mass of the current optical companion 
from ~ 2.9Mo to ~ 5.7 Mo, and the binary's probability 
of remaining bound after the supernova from 38% to 62%. 
With our adjusted prescription for the helium core mass of 
the stripped companion and fwr = 0.2, Cygnus X-1 has a 
596 probability to merge within 14 Gyr as a binary BH. 


We thus find that it is possible that a small fraction of 
HMXBs like Cygnus X-1 could form merging binary BHs, 
although this conclusion is sensitive to the treatment of 
mass transfer from MS donors in population synthesis mod- 
els and to the natal kick distribution of relatively low-mass 
BHs. If systems like Cygnus X-1 do become progenitors of 
gravitational-wave events, this would impact the predicted 
spin distribution of merging binary black holes (Kushnir et al. 
2016; Zaldarriaga et al. 2018; Belczynski et al. 2017; Fuller 
& Ma 2019; Bavera et al. 2020). Gravitational-wave obser- 
vations could ultimately address this possibility by resolving 
the spin distribution with more events (e.g., Farr et al. 2017). 
Meanwhile, wide, non-merging binary BHs could potentially 
be observable through microlensing (Eilbott et al. 2017). 


5. CYGNUS X-1: CAVEATS AND CONUNDRUMS 


We show that the current properties of the Cygnus X-1 sys- 
tem imply a reduction in Wolf-Rayet wind mass loss rates for 
exposed HeMS stars. These results depend on several key as- 
sumptions. 

We assumed that the optical companion did not experience 
Roche lobe overflow in its past. This assumption is consis- 
tent with the challenge of detaching from mass transfer once 
it commences given that the companion has roughly twice 
the mass of the BH, as explained in section 2.4. However, it 
is somewhat surprising that several HMXBs with well mea- 
sured properties - Cygnus X-1, LMC X-1 and M33 X-7 - 
share not only a high BH spin, but also a similar evolutionary 
state. Selection effects favour observing bright, long-lived 
systems, i.e., those with massive main-sequence donors that 
are close to Roche lobe filling (enabling more efficient ac- 
cretion). There may also be an evolutionary stalling point, 
increasing the number of systems in this phase. 

The latter scenario could indicate that the systems do 
manage to detach and resume mass transfer multiple times. 
While this would negate our wind mass loss rate conclusions, 
it would imply that much less angular momentum is car- 
ried away during non-conservative mass transfer onto a BH 
than we expected. Assuming non-conservative mass trans- 
fer (valid if accretion onto a BH is Eddington-limited) from 
a donor that is twice as massive as the accretor, the specific 
angular momentum of the material ejected from the binary in 
units of the binary's specific orbital angular momentum must 
be y < 0.85 in order to avoid a decrease in the size of the 
Roche lobe (see Eq. 2.1 for the change in orbital separation). 
For comparison, isotropic re-emission from the BH corre- 
sponds to y = 2. Conversely, if y = 2, the companion could 
still disengage from mass transfer if its radius shrinks faster 
than the size of the Roche lobe in response to mass loss. This 
would require the adiabatic logarithmic derivative of radius 
with respect to mass to exceed ¢ = dlog R/dlog M > 1.54, 
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which may be possible for stars in the late phase of their MS 
evolution that have already lost some mass. 

It is also possible that the primary is not fully stripped dur- 
ing the mass transfer episode, but retains about ~ 0.1 Mo 
of its hydrogen envelope (e.g. Yoon et al. 2010; Yoon 2017; 
Bersten et al. 2014; Gótberg et al. 2017, 2018; Laplace et al. 
2020). The changed surface abundance could lead to reduced 
mass-loss rates until the remaining hydrogen is completely 
removed. However, it is not clear whether retaining an en- 
velope of a fraction of a solar mass could be sufficient to 
prevent Wolf-Rayet-like winds. In any case, whether Wolf- 
Rayet wind mass loss rates must be lower than anticipated 
or whether stars that experience mass transfer in binaries are 
only partially stripped, the impact on binary evolution is sim- 
ilar: there is less mass loss than previously assumed. In fact, 
our Wolf-Rayet wind reduction factors can be broadly inter- 
preted as constraints on winds from stripped stars, whether 
they are naked helium stars or retain a small hydrogen-rich 
envelope. 

Naked helium cores can expand significantly in the last 
stages of their lives, potentially leading to another mass trans- 
fer episode from the BH progenitor late in the evolution. 
However, the degree of expansion is very mild for stars with 
initial masses > 20M, at near-solar metallicities (Yoon et al. 
2010; Hirai 2017; Laplace et al. 2020), so Cygnus X-1 is un- 
likely to have experienced such mass transfer. 

As a consequence of its mass accretion history, the sec- 
ondary may be over-luminous relative to single stars of the 
same total mass (Dray & Tout 2007, but see Hellings 1983, 
who concludes that they MS accretors quickly return to 
single-star models, and Braun & Langer 1995, who reach 
the opposite conclusion and find that accretors are under- 
luminous). Since we use single star evolutionary tracks to 
estimate the properties of the secondary star, this can affect 
our wind constraint from the no-RLOF condition. However, 
since we choose a stellar model that matches the observed ra- 
dius of the secondary at the present day, we do not anticipate 
the impact to be significant. 

Finally, as figure 3 shows, the level of reduction in the 
winds is sensitive to the assumed metallicity of Cygnus X- 
1. We now discuss this in more detail. 

Shimanskii et al. (2012) report that helium, carbon, oxy- 
gen, aluminium, sulfur and iron are overabundant by [X/H]- 
0.23-0.43 dex compared to the solar values (Anders & 
Grevesse 1989). Nitrogen, neon, and silicon have an even 
higher overabundance of [X/H]=0.69-0.94 dex. These val- 
ues appear robust against variations due to orbital motion and 
Roche-lobe filling factors, although some hydrogen and he- 
lium lines are sensitive to variations in the wind (Shimanskii 
et al. 2012). 

Previous accretion from the BH progenitor could signifi- 
cantly alter the chemical profile on the surface of the com- 


panion. For example, the detailed models of Qin et al. (2019) 
predict the observed enhancement of companion nitrogen 
abundances as a consequence of late main sequence mass 
transfer from the BH progenitor. In addition to direct ac- 
cretion, which is expected to enhance helium and nitrogen 
abundances, the deposited angular momentum can lead to a 
dramatic spin-up of the MS star (e.g. Packet 1981), although 
spin-up to near break-up frequencies may suppress subse- 
quent accretion. The surface could then be enhanced by he- 
lium and CNO-elements due to rotational mixing (Meynet 
& Maeder 2000; Heger & Langer 2000; Przybilla et al. 
2010). The rotational mixing might also make the star over- 
luminous compared to a non-rotating model (Langer 1992). 
The optical companion is observed to be tidally locked at 
present, with an inferred ratio of the rotational to orbital fre- 
quency of 1.05 + 0.10 Miller-Jones et al. (2021). Assuming 
a present-day rotational period of 5.6 days, the rotational fre- 
quency is a third of the Keplerian (break-up) frequency at the 
stellar equator). Alternatively, as discussed above and con- 
trary to the channel assumed in this work, the MS compan- 
ion may have been partially stripped by mass transfer onto 
the BH or its progenitor after the initial mass transfer phase 
from the primary, revealing deeper layers of the star. 
Although these mechanisms could be responsible for the 
overabundance of some of the elements, they have difficulty 
in explaining the overabundance of late stage burning ele- 
ments such as silicon and iron. This implies that these el- 
ements were primordially enhanced or were deposited from 
the progenitor of the BH in the final stages of its evolution. 
High primordial abundances imply that both stars in the bi- 
nary had a high metallicity at birth, which therefore requires a 
very strong reduction in the mass-loss rate (a factor of ~ 10) 
following the constraints described in section 3. On the other 
hand, a weak explosion induced by the collapse of the core 
can lead to an ejection of a small fraction of the outer part 
of the envelope at very low velocities. Because most of the 
envelope is assumed to fall back into the BH, the ejected ma- 
terial will be barely above the escape velocity, and could be 
efficiently accreted by the companion in a RLOF-like man- 
ner. If the heavy elements synthesized at the centre are effi- 
ciently mixed up to the outer regions before the inner slower 
material starts falling back, these elements can accrete onto 
the surface of the secondary. Such efficient mixing of heavy 
elements has been observed in supernovae such as SN1987A 
and Cassiopeia A (e.g. Utrobin et al. 1995; Fesen et al. 2006), 
and has been reproduced in 3D supernova explosion simula- 
tions (e.g. Hammer et al. 2010; Wongwathanarat et al. 2015, 
2017), while Liu et al. (2015); Hirai et al. (2018) explore 
the contamination of a MS companion by supernova ejecta. 
However, it is not clear whether similar degrees of mixing 
can be induced in failed supernovae that form BHs rather 
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than neutron stars, and the abundance pattern of the Cygnus 
X-1 companion merits further investigation. 

Regardless of whether we assume that the observed com- 
panion metallicity of Z = 0.03 (Shimanskii et al. 2012) is 
primordial, or use the lower metallicity of Z = 0.01 based 
on HD 227460 (Daflon et al. 2001), we conclude that the ob- 
served properties of Cygnus X-1 require a reduction in Wolf- 
Rayet winds to ~ 5—40% of their previously assumed values 
in the context of our assumed evolutionary channel. 

Recent theoretical modelling of mass-loss from stripped 
stars (Vink 2017; Sander et al. 2020; Sander & Vink 2020) 
points to reduced mass-loss rates compared to earlier litera- 
ture. Moreover, these models suggest a steep dependence on 
the Eddington factor, which can change significantly during 
the lifetime of the stripped star. This indicates that extrapo- 
lating the empirical Wolf-Rayet mass-loss rates to the entire 
duration of the stripped star life is misleading. Our results 
are qualitatively consistent with these findings. The reduced 
mass-loss rates could also be attributed to strong wind clump- 
ing, which is expected to occur in line-driven winds due to 
radiative instabilities (Owocki et al. 1988; Sundqvist et al. 
2018). Clumping of winds has been indirectly observed for 
massive MS stars in X-ray binaries (El Mellah et al. 2018; 
Lomaeva et al. 2020) and stripped stars may also experience 
high degrees of clumping. 


We further find that HMXBs like Cygnus X-1 form 
through a different evolutionary channel than the bulk of 
merging binary black holes (see, e.g., Mandel & Farmer 
2018, for a review). However, a fortuitous natal kick accom- 
panying the birth of the secondary BH could lead Cygnus X-1 
to merge as a BH binary within 14 Gyr. Gravitational-wave 
observations may be able to constrain the contribution of this 
channel to the formation of merging binary BHs through spin 
measurements. 
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APPENDIX 


A. WIND MASS-LOSS RATE FOR STRIPPED STARS 


In this appendix, we describe the particular parametrized formalism we used to model wind-driven mass loss from stripped 
stars. While more recent theoretical models are available (e.g., Vink 2017; Sander et al. 2020), this provides us with a convenient 
framework for investigating the empirical constraints placed by Cygnus X-1. 

Hamann et al. (1995) formulated a prescription for the wind mass loss rate for stripped or Wolf-Rayet stars as a function of the 
mass and luminosity of the stripped star. Hamann & Koesterke (1998) expanded on this work by reducing winds by a factor of 
VD where D is the wind clumping factor (Moffat et al. 1988; Nugis et al. 1998). Vink & de Koter (2005) further introduced a 
metallicity dependence to the winds. Combining the effects of clumping and metallicity leads to the prescription 


1. 
(dM /dt)wr = elt Mo yr, (Al) 


where term | is the result of Hamann et al. (1995); Hamann & Koesterke (1998) and term 2 is from Vink & de Koter (2005). 
Setting D = 100, i.e., reducing winds by a factor of 10, recovers Eq.(9) of Belczynski et al. (2010) and is consistent with the 
winds of Yoon et al. (2010). 

We follow Barrett et al. (2018) in scaling the prescription of Belczynski et al. (2010) by a multiplicative factor fwr in order to 
parametrise the uncertainty in the wind mass-loss rates: 


L 15 g %86 
BH A Mo wi. 


A2 
Io X (A2) 


(dM/dt)wr = fwr x 107 
Note that this is not the same fwr as used in Yoon et al. (2010) because ours already assumes a reduction of the original wind 
prescription of Hamann et al. (1995). The default assumption of fwr = 1 corresponds to the default models of Belczynski et al. 
(2010); Yoon et al. (2010); Stevenson et al. (2017); Qin et al. (2019). 

It is challenging to interpret our constraints on fwr in terms of a wind clumping factor D. A direct interpretation of fwr = 0.2 
would imply a clumping factor of D = 2500 in the model of Eq. (A1). Nugis et al. (1998); Nugis & Lamers (2000) used 
a clumping factor ranging from 10 to 30; D could be as high as 16 according to Hamann & Koesterke (1998). More recent 
theoretical models by Sander et al. (2017) use depth-dependent clumping factors and suggest a maximum value at infinity of 
Ds, = 10 to be consistent with observations of electron scattering wings. In later works however, Sander et al. (2020) and 
Sander & Vink (2020) propose 50 as a maximum upper limit for Dao, by comparison with previous theoretical works (Gráfener 
& Hamann 2005) and O/B star analyses (e.g. Bouret et al. 2012; Mahy et al. 2015). A clumping factor of 2500 seems extra- 
ordinary compared to previous models. As mentioned in the Section 5, additional constraints on the clumpiness of stellar winds 
from massive stars can be obtained from X-ray binaries (e.g. Lomaeva et al. 2020; El Mellah et al. 2018; Grinberg et al. 2017). 

In light of the above, our reduction of Jup is probably best interpreted as an overall constraint on mass loss from massive 
stripped stars (Vink 2017), perhaps indicating a different dependence on metallicity or luminosity (cf. Sander et al. 2020; Sander 
& Vink 2020), rather than a specific change in the clumping. 


